Reversible to Irreversible Flow Transition in Periodically Driven Vortices 
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We show that periodically driven superconducting vortices in the presence of quenched disor- 
der exhibit a transition from reversible to irreversible flow under increasing vortex density or cycle 
period. This type of behavior has recently been observed for periodically sheared colloidal suspen- 
sions and we demonstrate that driven vortex systems exhibit remarkably similar behavior. We also 
provide evidence that the onset of irreversible behavior is a dynamical phase transition. 
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Recent experiments on periodic shearing of colloidal 
suspensions in a viscous media where thermal effects are 
negligible have shown a transition from reversible behav- 
ior to an irreversible motion as a function of increasing 
particle density [J, The reversible and irreversible 
regimes were identified by measuring the net displace- 
ment of the colloids after each cycle of a periodic shear. 
In the reversible regime, the colloids return to their initial 
position at the end of each cycle, while in the irreversible 
regime, the colloids do not return to their starting point. 
A stroboscopic measurement of the displacements in the 
irreversible regime reveals that the colloids undergo an 
anisotropic random walk, with larger displacements in 
the shear direction. For a given colloid density, there is a 
threshold for the transition from reversible to irreversible 
behavior as a function of the strain amplitude or the 
distance the particles are sheared per cycle. The strain 
threshold decreases as the colloid concentration increases, 
indicating that colloid-colloid interactions play an impor- 
tant role in the transition. More recent modeling of this 
system has provided evidence that the reversible to ir- 
reversible transition is a nonequilibrium phase transition 
with power law divergences near the transition [I]. Other 
experiments on dilute sheared colloidal suspensions also 
indicate the importance of the particle interactions in 
producing irreversible or chaotic flow behaviors 

In this work, we consider whether the general features 
of the reversible and irreversible behaviors observed in 
the sheared colloidal system can be realized in a wider 
class of nonequilibrium many-particle systems. Driven 
particles moving over random quenched disorder provide 
an ideal class of systems for studying this issue [B| . Physi- 
cal realizations of such systems include vortices in type-II 
superconductors @, 01 , driven Wigner crystals Q , mag- 
netic bubble arrays [9(, driven pattern forming systems 
flpj ] . and colloids moving over random landscapes 11|. 
In the presence of strong quenched disorder, plastic flow 
regimes occur in these systems even in the absence of any 
thermal effects. In the case of vortices, it has been es- 
tablished that transitions from plastic flow to smectic or 
elastic flow states can be induced as a function of increas- 
ing dc driving force [f|. It is not known whether there 



could also be a reversible to irreversible flow transition 
in the presence of some form of periodic forcing. 

To address this question, we examine the specific sys- 
tem of vortices interacting with random quenched disor- 
der, and measure vortex displacements after successive 
cycles of a periodic drive. Vortex lattices have been used 
extensively as a general system in which to understand 
various nonequilibrium many-body effects since numer- 
ous simulations and experiments can be performed read- 
ily over a wide range of parameters. It has been shown 
that many of the nonequilibrium effects observed in the 
vortex system also occur in other systems such as driven 



colloidal assemblies [llj, [12[; thus, we expect that our 



results will be generic to other driven interacting point 
particle systems with quenched disorder. 

The simulation system consists of a two-dimensional 
sample of size L x L with periodic boundary conditions 
in the x and y directions. The sample contains N v vor- 
tices at a vortex density of n v = N v /L 2 . The vortex- 
vortex interaction force is F™ = Yl^j foKi(Rij / , 
where K\ is a modified Bessel function, iJy = |Rj — Rj| 
is the distance between vortex i and j located at R^ 
and Rj, = (R, - Rj)/Rij, f = o /(27r^ o A 3 ), 

A is the London penetration depth, and <p = h/2e is 
the flux quantum. Lengths are measured in units of 
A and the sample size is L = 64A. For Rij > A 
the vortex interaction falls off rapidly, so a finite cut- 
off on the interaction range is placed at Rij = 6A. The 
quenched disorder is modeled as N p — 100 randomly 

= 0.3A and 



placed parabolic pinning traps of radius r 
strength F p , with F? = £f p FP(R lk /r p )Q(r p - R, k )R lk , 
where 8 is the Heaviside step function, Ri k = |R; — R?|, 
Rifc = (Rj — Rj^/Ry-, an d -R-fe i s the location of pin k. 
The periodic drive is produced by the Lorentz force from 
a square wave current of period r, which we model as 
F ext (t) = F ex *sgn(sin(27rt/T))x, such that a vortex un- 
der the influence of only the driving force undergoes a 
displacement of d = F ext T/2rj in half a drive cycle. In 
this work, except where indicated, we fix F ext = 0.8/o 
and vary d by changing r. The overdamped equation of 
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motion for a single vortex i is 
dRi _„ _ 



F e:Et (t), 



(1) 



where 77 = 1 is the viscous damping term. The vortices 
are initialized in random positions. We use the same 
protocol as in Ref. [l[ to stroboscopically determine the 
vortex locations after each drive cycle and extract the 
effective diffusivities. The mean square vortex displace- 
ment R x and R~ in the x and y directions is measured 
as a function of the number n of drive cycles: 



N. 

R d a(nr) = N- 1 $}(Ri(t + nr) - R,(t )) • a} 2 (2) 



where a = x,y. From this quantity we follow Ref. [l[ and 
define an effective diffusivity D a — R~ t /2t. When the be- 
havior is reversible the vortices return to their initial po- 
sitions after each drive cycle and R x y and D x<y are zero. 
We observe an initial transient period of finite R%. y , af- 
ter which the system settles into a stationary state that 
is either reversible or irreversible. The diffusivities are 
measured only after this stationary state is reached. A 
similar transient behavior was observed in the shearing 
experiments of Ref. [l[ . We also measure the vortex ve- 
locity in the x direction given by V x — J^i " v, ■ x, 
where Vj is the velocity of vortex i. In all of our simula- 
tions, N v > N p so that we avoid any regimes where the 
vortices are all pinned. For all of the parameter sets con- 
sidered in this work, the vortices in the reversible regime 
are moving during the drive cycle. 

We find a transition from reversible to irreversible be- 
havior as a function of displacement per cycle d and vor- 
tex density n v . In Fig. 1(a) we highlight the locations 
of the vortices every O.Olr during n = 4 four cycles with 
d = 48 A and n v = 0.052/ A 2 . Most of the vortices are 
moving while a small number of vortices remain pinned. 
In Fig. 1(b) we illustrate the same system with the vortex 
positions marked once per cycle during n = 100 cycles. 
Here, all of the vortices return to their initial position 
at the end of each cycle, so all of the points marking 
the position of an individual vortex coincide. This in- 
dicates that the system is undergoing a reversible flow. 
For increasing d there is a transition to an irreversible 
flow state. In Fig. 1(c) we show the stroboscopic loca- 
tions of the vortices in n = 25 cycles for the same system 
as in Fig. l(a,b) but with d = 160A. In this case, the 
vortices do not return to their initial positions after each 
cycle. When the number of cycles increases, the vortices 
continue to move further from their starting point after 
each cycle. Figure 1(d) illustrates the anisotropic nature 
of the vortex displacements for the system in Fig. 1(c) 
during n = 100 cycles. 

In Fig. 2(a) we plot the cumulative mean squared dis- 
placements R x and R y versus drive cycle number for the 
sample in the irreversible regime shown in Fig. l(c,d). 
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FIG. 1: Black dots: Stroboscopic vortex positions marked 
every O.Olr in (a) and every r in (b,c,d). (a) Vortex positions 
during n = 4 cycles in a sample with n v — 0.052/A 2 and 
d = 48A. (b) Vortex positions in the same sample as (a). In 
this reversible regime, all of the vortices return to their initial 
positions after each cycle, (c) Positions in a sample with a 
longer d — 160A in n — 25 cycles. The vortices do not return 
to their initial positions and the behavior is irreversible, (d) 
The same sample as (c) for n = 100 cycles showing that the 
displacements are anisotropic. 



The dashed line fits indicate that the vortices are per- 
forming an anisotropic random walk similar to that found 
in the irreversible regime for the sheared colloidal system 
of Ref. [l[ . By conducting a series of simulations in which 
we vary d and n v , we can extract the effective diffusivities 
and determine the thresholds for the irreversible behav- 
ior. 

In Fig. 3(a) we plot D x and D y versus d for a sam- 
ple with fixed n v = 0.052/A 2 . There is a clear threshold 
d c = 80A marking the transition from reversible behavior 
with D x — D y — to irreversible behavior with nonzero 
D x and D y . Above d c , D x and D y monotonically in- 
crease with d. Figure 3(b) shows D XjV for fixed d = 96A 
and varied vortex density n v . There is a threshold den- 
sity n v — 0.034/A 2 below which the system is reversible 
and above which the behavior becomes irreversible. The 
diffusivity passes through a peak as n v increases above 
the threshold, while D x and D y decrease at higher values 
of n v when the increasing strength of the vortex- vortex 
interactions leads to an effective caging effect. The inset 
of Fig. 3(d) illustrates the appearance of two additional 
phases at higher n v . For 0.3/A 2 < n v < 0.36/A 2 , D y = 
but D x remains finite, resulting in a smectic phase Q 
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FIG. 2: (a) R x (top line) and Ry (bottom line) versus the 
number of drive cycles n for the system in Fig. l(c,d) in the 
irreversible regime. The particles undergo an anisotropic ran- 
dom walk motion. Dotted lines: fits used to obtain D x and 
D y . (b) S(f) for the irreversible regime in Fig. 1(c) (top line) 
and the reversible regime in Fig. 1(a) (bottom line). 
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FIG. 3: Diffusivities D x and D y versus d for fixed n v — 
0.052/A 2 . (b) D x and D y versus n v for fixed d = 96A. (c) D x 
and D v versus F p for fixed n v = 0.052/A 2 and d = 160A. (d) 
Threshold displacements d c for the transition to irreversible 
behavior versus n v for fixed F p — 1.25/o. Inset: D x (circles) 
and D y (squares) at n v > 0.25/A 2 continued from (b). 



where the vortices move in one-dimensional chains that 
are decoupled in the y direction. For n v > 0.36/A 2 , the 
system enters a jammed phase where the vortex- vortex 
interactions dominate, the vortices form a rigid lattice, 
and the motion is reversible again. 

In Fig. 3(c) we illustrate the effect of the pinning 
strength F p on the dynamics for a sample with d = 96A 
and n v ~ 0.052/A 2 . For weak F p the system acts re- 
versibly, since when F ext > F p , all the vortices depin 
and there is no local shear in the vortex lattice. In the 
colloidal shearing experiments of Ref. |l( , the shear ge- 
ometry generated a velocity differential in the system, 
and when the total shear was too small, the colloids 
moved reversibly. At high F p , the vortices are immo- 
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FIG. 4: (a) V x versus cycle number n for (light line) the 
reversible regime in Fig. 1(a) and (heavy line) the irreversible 
regime in Fig. 1(c). (b) D x from Fig. 3(a) vs (d — d c )/d c . 
Dashed line: logarithmic fit D x oc (d — dc) 13 with /3 = 0.59. 
(c-f) Overlapping time traces of V x during 10 cycles plotted 
versus t' = t — t mod r. Reversible regime from Fig. 1(a): 
(c) lower half cycle, (d) upper half cycle. Irreversible regime 
from Fig. 1(c): (e) lower half cycle, (f) upper half cycle. 



bilized by the pinning sites and the behavior is reversible 
again. Irreversible behavior occurs for intermediate val- 
ues of F p , as shown in Fig. 3(c) by the peak in D x and 
D y at F p w 1.5/ . In Fig. 3(d) we plot d c , the thresh- 
old displacement per cycle required to induce irreversible 
behavior, versus vortex density n v , showing that d c de- 
creases with increasing n v . There is a divergence in d c 
near the matching density n v — 0.0244/A 2 where the 
number of vortices equals the number of pinning sites. 



If N v < N p , all the vortices are pinned and the sys- 
tem enters a reversible state. As N v — » N p from above, 
the distance between nonpinned vortices diverges and the 
vortices must be driven over longer distances in order to 
permit the nonpinned vortices to interact with each other 
and produce the irreversible behavior. 

In general, our results are in good agreement with the 
shear experiments of Ref. [l|, where a threshold strain 
which decreases with increasing colloid density must be 
applied to produce irreversible motion. d c diverges as the 
threshold is approached from above. Our results suggest 
that if the colloidal experiment were performed at higher 
colloid density than considered in Ref. [lj , the diffusivity 
should decrease dramatically and could drop to zero if a 
jamming transition occurs. It may also be possible that 
the colloidal motion will become banded, which would be 
analogous to the smectic state. 

In the shear experiments the transition from reversible 
to irreversible behavior can be visualized directly by 
video microscopy while for most vortex systems, di- 
rect imaging of the vortex motion is difficult. A straight- 
forward method for probing the nature of the vortex dy- 
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namics is through voltage noise fluctuations, where the 
voltage is proportional to the vortex velocity. In Fig. 4(a) 
we plot V x versus cycle number n for the same system 
shown in Fig. l(a,b) in both the reversible and irre- 
versible regimes. The velocities are sampled at a rate 
of O.Olr. In Fig. 4(a), we see that V x follows the square 
wave of the driving force and that the average V x for 
the two states is similar. We can superimpose the veloc- 
ity curves in each cycle by plotting V x versus t', where 
if = t — t mod t, and compare the velocity fluctuations 
in the two regimes. In Fig. 4(c,d) we plot the super- 
imposed fluctuations of V x for 10 cycles in the reversible 
regime for the negative and positive portions of the cycle. 
Here, the fluctuations of V x are nearly identical in each 
cycle, indicating that the vortices are flowing along the 
same trajectories during each cycle. In the irreversible 
regime, shown in Fig. 4(e,f), the velocity fluctuations dif- 
fer during each cycle and the curves do not overlap. We 
observe a similar behavior for the fluctuations in the y- 
direction. Since the power spectrum of the fluctuations 
can be measured readily in experiment [IH , we compute 
S(f) = \j\V x {t)\e- 2 ^ t dt\ 2 . Fig. 2(b) shows that a 
broad spectrum appears in the irreversible regime, while 
in the reversible regime the noise power is much lower 
and strong harmonic frequencies appear. These results 
show that the reversible to irreversible transition can be 
identified through noise measurements. A more direct 
visualization experiment could be performed by periodi- 
cally driving a two-dimensional colloid system interacting 
with a random substrate [12] or optical traps JtI- 

The colloidal shearing experiments of Refs. [l|,|3| found 
evidence of a dynamical phase transition into the irre- 
versible state, along with an initial transient behavior of 
a duration that diverges near the transition. We find a 
similar transient behavior in our system. If we consider 
D x and D y in Fig. 2(a) as an order parameter, Fig. 4(b) 
shows that we find a reasonable scaling of D x oc (d — d c )^ 
with [3 = 0.59, consistent with either directed percolation 
in two dimensions [IH or conserved directed percolation. 
The latter was recently proposed as the universality class 
of the colloidal shear experiments [lj| . It would be inter- 
esting to determine whether a similar power law diver- 
gence occurs in the average displacements near the strain 
threshold in the colloidal system. 

In summary, we have shown with numerical simula- 
tions that a non-thermal transition from reversible to 
irreversible flow behavior occurs for periodically driven 
vortex systems in the presence of quenched disorder. The 
transition is characterized by performing stroboscopic 
measurements of the vortex locations after each drive cy- 
cle. We find that there is a threshold in the displacement 
per cycle above which the system acts irreversibly. In the 
irreversible regime, the vortices undergo an anisotropic 
random walk. The threshold decreases as the vortex den- 
sity increases. The transition can also be identified by 
analyzing the voltage noise fluctuations during each cy- 



cle. Our results are remarkably similar to recent experi- 
mental observations of a transition to irreversible flow for 
sheared colloids. In the colloid system, the velocity dis- 
persion created by the shearing process permits neighbor- 
ing colloids to undergo different random displacements. 
In the vortex system, the pinning sites create velocity 
dispersion by generating local shear in the vortex lattice. 
The transition to the irreversible state is consistent with 
either two-dimensional directed percolation or conserved 
directed percolation. Our results suggest that the behav- 
ior of the colloidal shear experiments may be general to 
driven particle systems with quenched disorder. 
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